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We construct the dilaton potential in the gravity dual theory of AdS/ QCD 
for the warp factor of Refs. J7], Using this AdS^-metric with properties 
similar to QCD, we find that the gravity dual leads to a meaningful gauge 
coupling in the region between the charmonium and bottonium mass, but 
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• differs slightly from QCD in the extreme UV. When we fix the ultraviolet 

behavior in accord with the ^-function, we can obtain good agreement with 
the overall heavy quark- antiquark potential. Although the leading order 
^ ! proportional to — differs from perturbative QCD, the full potential 

CN . agrees quite well with the short distance QCD potential in NNLO. 
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1 Introduction 

o\ . 

Maldacena's conjecture [5] states that, in the low energy limit, the large- N c Af = 4 super 
Yang-Mills field theory in four-dimensional space is equivalent to the type I IB string theory 
in AdS§ x S 5 - space. This may provide a possibility to solve the longstanding problem of 
strongly coupled QCD in the low energy limit. But a rigorous top-down approach is still 
far from giving any experimental prediction. Using the idea of holography, bottom-up 
models can quantitatively reproduce many results of QCD in the low energy limit gained 
by other existing methods like lattice QCD. The core of these bottom- up models is to find 
a reasonable non-conformal metric of the AdS5-spa.ce, which incorporates relevant physical 
information. 

To obtain such a metric, one can either maintain the basic form of the conformal AdS$- 
metric and introduce a factor by hand, or assume a very general form of the metric, and 
calculate its explicit parameters. In this paper, we further investigate the warp factor 
given in Refs. pQ [2] based on the interpretation of the fifth-dimension z as a coordinate 
proportional to the inverse energy resolution. The warp factor of Refs. [TJ [2] resembles the 
running coupling in QCD - having a strong inverse logarithmic growth in the infrared. The 
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1 Introduction 



Minkowskian form of the metric with the warp factor is 

ds 1 = h(z)—^—r(-dt 2 + dx 2 + dz 2 ), (1) 

with 

h{z) = (2) 

This metric with asymptotically conformal symmetry in the UV and infrared slavery in 
the IR region yields a good fit to the heavy QQ-potential with 

A = 264 MeV, (3) 
e = A 2 l 2 = 0.48. (4) 

In Ref. [2], we have used this metric to calculate the expectation value of one circular 
Wilson loop (W) and the correlator of two concentric circular Wilson loops (WW) from 
the Nambu-Goto action of the form: 



Sng = J d 2 ^det hob, (5) 

where l s is the string length defined above. The induced world-sheet metric in the Nambu- 
Goto action is called 7t &: 

dX»dX» 

Hab ~ ^ a ^ . [b) 

This calculation gives reasonable results for one and two Wilson loops. In addition to 
the confinement physics, we have derived the gluon condensate at zero temperature. The 
above metric, Eq. flTJ, is given in the string frame, and all the calculations mentioned above 
are done in this frame. Low energy string theory can be approximated further by a gravity 
action with a background scalar field, the dilaton <fi. When the conformal symmetry is 
broken, the gravity action contains a dilaton potential V(4>), which is no longer constant 
and can be calculated from the metric with the help of the Einstein equations. Relating 
the energy scale of the gauge theory living on the boundary of the AdS 5 -spa.ce and the 
bulk ^-coordinate, we can predict the running gauge coupling inside the region between 
the charmonium and the bottonium mass. 

When we check whether this solution of the running coupling is compatible with the 
behavior of the QCD coupling in the far ultraviolet, we find small deviations. The reason is 
that our metric in spite of many good phenomenological features has conformal invariance 
in the ultraviolet. By modifying the dilaton potential V(<f>), we can make the potential 
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consistent with the QCD /3-function and the heavy QQ -potential. We find good agree- 
ment with the short distance behavior of the QQ-potential. Finally we close the circle of 
investigation and recalculate the warp factor h(z) of the modified metric in string frame. 

The outline of the paper is as follows: After the introduction in Section [TJ Section [2] 
describes the formal construction of the dilaton potential from the metric in the string 
frame. Section [3] shows how the energy scale of QCD is related to the bulk coordinate z 
and gives the numerical calculation of the dilaton potential. In Section HJ the resulting 
QCD running coupling is derived and the ultraviolet behavior of the potential is improved. 
Section [5] applies the UV-improved dilaton potential to a calculation of the heavy quark 
potential and the glueball spectrum. We give in Section [6] a discussion about the modified 
dilaton potential, the running coupling and the new warp factor. Finally, Section [7] gives a 
final discussion and our conclusions. In Appendix A the infrared and ultraviolet properties 
of the 5-dim Nambu-Goto theory with the metric of Eqs. (JT|) and (T5]) is investigated. In 
Appendix B we derive the infrared properties of the UV-improved dilaton potential. In 
Appendix C we show technical details for the analytical computation of the heavy QQ- 
potential in the small distance regime. 

2 Construction of the Gravity Dual Theory with Dilaton Potential 

We argued in the previous section that the warp factor given by Eq. flTJ proposed in Ref. [1] 
produces reasonable results, as shown in Refs. [TJ [2]. But the metric has to be consistent 
with gravity. In the low energy limit, string theory can be approximated by its gravity 
dual theory. In a top-down approach one obtains the II Aj II B effective action for the 
long-range fields of _D 3 -branes: 



Here R is the Ricci-scalar of the gravitational field, is the dilaton field and F 5 is the 
5-form flux originating from the branes. This string frame action is characterized by the 
exponential dilaton dependence in front of the curvature scalar. In order to keep the theory 
simple, we neglect the axion field a and other eventual space filling branes related to quark 
dynamics. After inserting the equation of motion for F§ back into the action, we transform 
the action from the string frame to the Einstein frame, in which the Einstein term has the 
conventional form by a Weyl rescaling of the metric: 



SlOD-Gravity 




d 10 xV^& e" 2 * (R + 4(V0) 2 ) 




(7) 



G E (X) = e-l*G*(X) 



(8) 
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2 Construction of the Gravity Dual Theory with Dilaton Potential 



Integrating over the Ss-space and combining several terms into a dilaton potential V(4>), 
we obtain the following five- dimensional action in the Einstein frame [I]: 

S 5 D~GravUy = ^| J d 5 Xy^G^ (r - ^d^<j> ~ V{<j>^j . (9) 

In principle, there can be an additional Gibbons-Hawking term in the action [51 [6] . But 
since this term does not contribute to the variation with respect to the metric G^ u , it 
does not affect the Einstein equations, which determine all the physical quantities at zero 
temperature. We must emphasize that, at finite temperature, only Einstein equations 
themselves are not enough to determine all the physical quantities, but the action is also 
important, hence at finite temperature the Gibbons- Hawking term contributes to the ther- 
modynamic quantities. In this paper, we only focus on the T = case, therefore, we can 
neglect the Gibbons-Hawking term. 

We assume that the background dilaton potential denoted by V(<p) incorporates some 
information of the 5-form F§ and higher curvature corrections extending the range of 
applicability. The background field <fi(z) has a z-dependence reflecting necessary corrections 
at higher energies. Let us comment on the dimensionalities of the quantities introduced 

in the action, Eq. ([9]). The Ricci scalar R has dimension \R] = — ^-ttt, k% has dimension 

length 

[ac§] = length 3 . Consequently, the bulk field <p is dimensionless, and the dimension of the 
dilaton potential V((j>) is [V(0)] = j — ^h 7 ' ^ ur s P ace_ time metric has Minkowski signature 
with the sign convention (— , +, +, +, +). We emphasize that all quantities in Eq. (Q have 
to be in the Einstein frame. In the remaining part of this paper, we simply write G^, 
instead of Gz v . After variation of the Einstein-frame action, Eq. (Q, with respect to the 
metric G^ v , one obtains the following equations of motion 

(R, v - \rg, v ) - (\d,<pd v <p - ^G^d^d^ + V{<f>))) = . (10) 
* i ' ^ — 2- ' 

The Einstein equations contain the energy momentum tensor T^, 

E^ u = T^ u . (11) 

With the above chosen signature of the metric and assuming just a flat metric, we obtain 
^oo — |(0') 2 + j where the prime denotes derivative with respect to z. Up to a 

normalization factor of the kinetic term, this agrees with the definition of the energy of a 
scalar field. 



1 Note in Refs. [HE] V{(f>) is defined as the negative of our potential. 
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The warped metric in the string frame, Eq. (OQ), becomes the Einstein-frame metric 
through Eq. QgJ): 

ds 2 E = e-^h{z)j^-^{-dt 2 + dx 2 + dz 1 ). (12) 

The resulting Einstein equations, Eqs. ffTUj) . can be solved for the dilaton. Besides the 
Einstein equations, the action given by Eq. (J5J) yields another Euler-Lagrange equation, 
which can be obtained by varying the action with respect to <p: 

This equation contains only redundant information and is not independent of the Einstein 
equations, Eqs. ( jiop . Hence Eq. (fT5|) needs no further treatment. This is a consequence of 
the fact that is actually not an independent field, but part of G^ v . 

Now we turn to the Einstein equations, Eqs. ([TO]) . To be consistent with Refs. [U [7], 
we use the following form of the space-time metric in the Einstein frame 

G^ = e 2 ^-diag(-l,l,l,l,l). (14) 

Comparing this metric with Eq. (|T2|) . we have 



,2A{z) 



e-HM^- (16) 

We first calculate the tensor E^ u in terms of A(z) and its derivatives. Then we compute 
the components of the energy-momentum tensor T^ v under the assumption that <p — 

- T u = T 22 = T 33 = T 44 = ~{<P>f - \e 2A ^V{cj>) : (16) 

T 55 = 2 -(<P') 2 -le 2A ^V(<P)- (17) 
Thus, we end up with only two independent equations of motion, namely 

3((A'(z)) 2 + A"{z)) = - 2 -(<f>T - \e 2A{z) V{<P\ (18) 

6(A'(z)) 2 = ^'f - \e 2A ^V{4>). (19) 
Adding these equations we obtain a formal expression for the dilaton potential: 



V(<f)(z)) = -e- 2A{z \9(A'(z)) 2 + 3A"(z)). 



(20) 
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3 The Energy Scale E(z) and the Solution for the Dilaton Potential 



Multiplying Eq. ( TT8|) by (—1) and then adding it to Eq. (|T9|) we find an important relation 
between the dilaton and the metric profile: 

{<P'? = \{{A{z)y-A\z)). (21) 

These two equations, Eq. ( 1201) and Eq. ( 12T1) . are identical with Eq. (2.16) in Ref. pi], up 
to a different sign of the dilaton potential. The different sign stems from the fact that 
we have a minus sign in front of the potential in the action, Eq. (Q. We have chosen the 
potential in this way, because we want the Too-component of the energy-momentum tensor 
to be the sum of the kinetic and the potential energy. We see that Eq. (T2~T]) depends on 
the profile A(z), which is a function of the warp factor h(z) and the dilaton field <fi(z) (cf. 
Eq. ffT5]) ). The resulting second order differential equation for <ft(z) needs two boundary 
conditions, which we will obtain from the QCD running coupling constant once the bulk 
coordinate z is connected with the energy scale E(z). 



3 The Energy Scale E(z) and the Solution for the Dilaton Potential 

Before we turn to the numerical solution of Eq. (12T1) . we must do some preparations. Per- 
turbative string theory is based on a topological expansion in the string coupling g s , which 
generates factors e~ x ^ in the string partition function, where x is the Euler characteristic 
of the string surface, e.g. \ = 2 for a sphere. Therefore the string coupling is proportional 
to e^. We choose the proportionality constant equal to unity, i.e. 

g s = e*. (22) 

The elastic amplitude of two closed strings (glueballs) of order 0{gf) corresponds to a 
contribution 0(a 2 ) in the boundary field theory. In general, the AdS/CFT correspondence 
relates the string coupling g s of perturbative string theory to the Yang-Mills coupling a in 
the following manner: 

9yM /no\ 

g s = — =a. (23) 

In order to simplify the notation, in the following we will use a to denote both the Yang- 
Mills coupling related to the string coupling and the QCD running coupling. 

The warping of the bulk space relates the bulk coordinate z to the energy scale E(z) 
associated with z via the gravitational blue-shift. It is pedagogical to use the radial co- 
ordinate r oc 1/z. In SU(N), one has N _D 3 -branes as gravitational source located at 
E r ^oo denotes the value of the energy scale at the holographic boundary 



In our notation the radius of the AdS5-spa.ce is L = ^, see Eq. (UJ). 
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r — > oo, when E r is the value of the energy at an arbitrary value of r. The blue-shift is 
given by the dimensionless ratio 



E r IGttir -> oo) 



where G 4t denotes the temporal component of the metric. In the limit r — > oo, we are far 
away from the branes, where the space-time is asymptotically flat, which yields G#(r — > 
oo) = —1. Hence, the blue-shift reads 

Er^oo = E r y/-G tt (r) or equivalently E r ^ = E z y/-G u (z). (25) 

In the unmodified AdS^-space, G u = — t^w and hence -E r _>oo = E z -^. For z — > 
0, we find the UV regime of the boundary field theory. This is in agreement with the 
intuition underlying the renormalization group interpretation of the ^-coordinate, which 
was instrumental to guess the warp factor h(z). For E z one can choose an arbitrary value 
of the energy scale. In order to simplify the expression, we choose the confinement scale 
given by A = 264 MeV. This leads to the following explicit formulas in the Einstein frame 



E r ^ = e-^^M_ (26) 



z 

= a-§^53 (27) 

z 

= e A{z) ■ A. (28) 

Suppose that we know the value of the coupling constant a at a given energy scale E = 
Er^oo, then we can find the corresponding value of z from Eq. ( ]2"7|h 

At a given value of z, 4>(z) = log (a) gives just one boundary condition for Eq. (1211) . In 
order to obtain the second boundary condition to Eq. (12~TT) . we need a second value of the 
QCD running coupling. It is not easy to choose two appropriate energy values. The reason 
is the following. In the UV limit it is questionable whether the AdS/CFT correspondence 
is still valid. On the other hand, in the IR limit, there is no reliable measurement of 
the running coupling. We think that the region between the charmonium mass and the 
bottonium mass, i.e. between 3 GeV and 8 GeV, is a reasonable region where the modeling 
of AdS/QCD with the warp factor of Eq. ([T|) should work well. Therefore, we propose as 
input values for a [8]: 



a(3 GeV) = 0.25241, a(8 GeV) = 0.18575. 



(29) 
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3 The Energy Scale E(z) and the Solution for the Dilaton Potential 



0(z) 



-3 



-2 



-1 



1 




z [GeV -1 ] 



Fig. 1 : The dilaton field profile 4> as a function of the bulk coordinate z, computed from the warp 
factor h{z) of Eq. 

This also means that we have implicitly set N c = 3 and Nf = 4. Although there is no entry 
in our model for color and flavor, we have fitted our final result to the Cornell potential, 
therefore we have implicitly made a choice for these parameters. 

The fact that one has to fix two integration constants using initial conditions contrasts 
with the analysis of Refs. [TJ [9], where the authors show that just one initial condition is 
enough. One possibility studied in Refs [TJ [9] is fixing one of the parameters by requiring 
that the bulk singularity is not of the "bad kind" , which means that the singularity should 
be repulsive to physical fluctuations. In our case the singularity in the infrared is of the 
"good kind". As we explain in Appendix A, the reason why we have to handle with two 
integration constants is that we do not use the perturbative /3-function as starting point. 

With the above conditions we solve Eq. (12T1) numerically, and obtain the dilaton field 
as a function of z (cf. Fig. [1]). Note the strong variations of the dilaton field for small 
and large z. For large z, infrared confinement at low energies is felt. 

In Appendix A we have investigated the infrared and ultraviolet properties of <p(z) 
analytically. Defining 



£ = zir - 



z . 



(30) 



with 



(31) 



we find in the infrared 




(32) 



9 



with 

co m = 4.55 GeV -1 , « m = -0.758 , (33) 

and in the ultraviolet 

(j){ z ) = _^HI + Kuv + (z) , ^^O. (34) 

The parameters k uv and related, in the sense that setting /% v in the UV then 

sets Km in the IR. u vv and u m are related in the same way. From a numerical computation 
of (j)(z) in the full regime < z <C z IR we find that Eq. (1331) , or equivalently Eq. 0291) , 
leads to: 

u vv = 0.1285 GeV" 1 , /% v = -1.386 . (35) 

We show in Fig. [2] the relation between the parameters k uv and k ir when co m and w uv are 
fixed to the values quoted in Eqs. ( 1331) and ( 1351) . The functional dependence is 

Kuv = Km - 0.628 . (36) 




Fig. 2: T/ie UV parameter k uv as a function of the IR parameter k [R , cf. Eqs. (32]) and \34\ ), 
when uj ir and oj uv are fixed to the values quoted in Eqs. \33\) and i35\) . The point corresponds 
K m = —0.758 and k uv = —1.386, which follows from the initial condition Eq. 129\) . The functional 
dependence is given by Eq. A36\) . 

Inserting <f>(z) into Eq. ( |26i) . we calculate immediately the associated energy scale 
E r ^oo{ z ), which is shown in Fig. [3] as a function of the z-coordinate in the fifth dimension. 
At the lower edge of the z-scale, we have a strongly increasing energy E^^z) (cf. Fig. [3]). 
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3 The Energy Scale E(z) and the Solution for the Dilaton Potential 



>oo(z) [GeV] 
60 




0.0 0.5 1.0 1.5 2.0 2.5 



z [GeV -1 ] 



Fig. 3: The energy scale as a function of the bulk coordinate z, corresponding to the warp factor 
Hz) ofEq. 



Consequently, the A{z) in the metric can be calculated from Eq. (128|) . With A(z), the 
other Einstein equation, Eq. (1201 . gives us V(z). Combining V(z) with </>(z), we obtain 
the dilaton potential V((f>), which is shown in Fig. HJ The result shows an approximately 



V(0) [GeV ] 




Fig. 4: The dilaton potential as a function of the dilaton field (p. The result corresponding to the 
warp factor h(z) of Eq. (Tj|) is displayed as a full (red) line. We show for comparison as a dashed 
(blue) line the modified dilaton potential given by Eq. (5$ with the parameters provided in Sec. [3J 

constant dilaton potential until <j>{z) = 1. Beyond this point the dilaton potential falls 
rapidly. Recall that in the AdS-sp&ce, the "cosmological" term is negative and slowly 
varying due to the asymptotic conformal behavior of the warp factor. In the conformal 
limit, 0' = 0, and the dilaton potential should have the value —j?, which can solve the 
Einstein equations, Eqs. f[2TJ|) and f[2"T|) . 

Questions about the stability of the vacuum because of the large and negatively un- 
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bound dilaton potential have to be analyzed, but it is well known that due to Breitenlohner- 
Freedman bound negative second order derivatives in the dilaton potential do not cause 
problems in the presence of gravity [TQl [TT] . 

4 Constraining the Dilaton Potential by the QCD -Function 

From the energy scale E(z) and the dilaton profile <p(z), we are now able to calculate the 
value of the strong coupling constant at any energy scale a = e^ z >. The gravity dual 
of string theory allows to interpolate the QCD coupling between our boundary values at 
3 GeV and 8 GeV in a satisfactory manner, as one sees from the comparison of a = in 
Fig. [5] with the strong coupling from the PDG web tool [Sj. The good description of the 
coupling adds another positive feature to the warp factor proposed in Refs. [H [2] . 




E [GeV] 



Fig. 5: The running coupling as a function of the energy scale. Full (red) line corresponds to the 
string theory result that follows from the relation a = e$ , using the warp factor h(z) of Eq. (TJ|). 
We show as a dashed (black) line the experimental values of the running coupling from the PDG 
data The two points correspond to the input conditions, Eq. $29\). 



When we investigate the behavior of the running coupling in the deep UV more closely, 
we expect some contradiction with QCD. As much as conformal behavior is favored in 
correlation functions, where the leading behavior up to logarithmic corrections is correctly 
reproduced, we have to deviate from the correct running of the coupling in the deep 
ultraviolet, since our metric h(z) has scale independence in this limit. We start with 
the definition of the /3-function: 



(37) 
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4 Constraining the Dilaton Potential by the QCD /3-Function 



Eq. (|28|) tells us that the energy can be expressed as the product of e A ■ A. Using a = e^, 
we obtain: 

da A d{e*) e+d<j> e^> ■ fjz) 
P = E dE =e A -d(^A) = ^A = A'(z) ■ (38) 

All quantities in the last expression are calculable from the warp factor h(z). In Fig. [6] 
we show the /3-function from our AdS/QCD model together with the QCD /3-function at 
two-loop level. In QCD the /3-function has the following form |12j : 



£(a) = -b a 2 - M 3 , (39) 



with 



b = — (—N c - -N f ) , and h = ^-( -N 2 - f-N c - —) N f ) . (40) 
2tt V 3 3 V 8tt 2 V 3 6 \ 3 NJ i ) v ' 

As argued before, we have set N c = 3 and iV/ = 4. In this case, 6o — |f , and b\ = j^z- 
One sees that the agreement is very good near a = 0.25, i.e. near the charmonium mass 
region, where the phenomenological adjustment has been done. For smaller values of a 
there is a sizable discrepancy. Explicitly in this region, the /3-function stemming from the 
warp factor oc —a, not oc —a 2 as in QCD. This explains the deviation visible in Fig. E] 
for small a. We have derived in Appendix A an analytic expression for f3(a) in the deep 
ultraviolet, cf. Eq. f ll 19j) . The strongly increasing warp factor h(z) in the infrared leads to 
a stronger beta function for large a. 

How can one repair this problem and make the dilaton potential consistent with QCD 
in the ultraviolet? The basic idea of Refs. [1] and [7] is to use the QCD /3-function itself 
as a starting point, and derive the metric from the /3- function. The resulting metric is 
then unambiguously consistent with the QCD /3-function, as expected, and the running 
coupling calculated from this new metric is necessarily correct. This procedure presents a 
systematic approach to define the dilaton potential in the ultraviolet. In the infrared, for 
large positive values of the dilaton field there remains the problem to choose a parametriza- 
tion of the potential. Calculations of the string tension have been proposed as tests of this 
parametrization at zero temperature [13] or of the spatial string tension at finite tempera- 
ture [H] . In our case, we will build on the phenomenological work done in Refs. (H [2] and 
will fit the constrained form of the potential to the heavy quark potential. 

In the following, we will review the important formulas given in Refs. jl] and [7]. In 
the so-called "domain wall coordinates" du = e A dz: 

ds 2 = e 2A (-dt 2 + dx 2 ) + du 2 , (41) 
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Fig. 6: The /3-function as a function of the running coupling a. We show as a full (red) line 
the string theory result that follows from the warp factor h(z) of Eq. {Hp, using Eq. (Tffffj) . We 
show for comparison as a full (blue) line the result corresponding to Eq. [51]) with the parameters 
provided in Sec. Dashed (black) line is the QCD (3-function given by Eq. A39\) . 



the Einstein equations become 

3A + 12A 2 = V, (42) 

A = ~4>\ (43) 

where the dot denotes the derivative with respect to u. After defining two auxiliary vari- 
ables 

W = ~A, (44) 

we may rewrite the Einstein equations using W: 

x 4 • dW Tr 4 fdW\ 2 64 TT r2 . . 

W{4>) plays the role of a superpotential. Several equivalent expressions hold for X: 

X = (47) 
da 



14 4 Constraining the Dilaton Potential by the QCD /3-Function 



When the /3- function is known, we can calculate X through the first expression, Eq. (1471) . 
With this X, W is obtained as a solution of Eq. (j45|) . Consequently, A and can be 
calculated through Eq. ( 1441) and Eq. ( 1481) . respectively. Finally, the general form of the 
dilaton potential V, is determined from the last equation of the three Eqs. (1461) as 

y (0) = V ■ (1 - X 2 ) • e-t (50) 

Therefore, for a /3-function given over the whole range of a, the dilaton potential is fixed. 
In the IR region, we do not know the correct form of the /3-function. If we want to impose 
confinement in the IR region, some forms of the /3-function are excluded, but there still can 
be several possible classes to achieve confinement [TJ. It is well known that the definition 
of the /3-function becomes dependent of the quantity one studies, when the gauge coupling 
becomes strong. Therefore, the question arises, which parametrization of the /3-function 
should one choose in the infrared. One possible choice is given by Eq. (5.1) in Ref. jl]. 
Here we propose another possible choice which combines the correct UV-behavior with 
some integrable form, 

(51) 

This parametrization has the required /3- function in the UV region as limit, 

a — > : /3(a) rs —boa 2 — b\a 3 , 

and confinement property in the IR region [7J . 

We can get a constraint on the parameters b 2 and a by demanding a good behavior for 
the running coupling. For a < a the QCD-coupling is strictly perturbative, whereas for 
a > a the /3-function is characterized by the non-perturbative linear term — b 2 a. When 
we consider the coupling in the region 0.6 GeV < E < 15 GeV, we obtain a good fit of the 
running coupling for values 1.2 < & 2 < 3 and a > 0.27 along the line § 

— = 5.09. (52) 
a 

The x 2 /d-o.f. is very close to its minimum in the entire region of 62- The perturbative 
running for energies larger than the charmonium mass is guaranteed by the limit on a. On 
the other hand, the /3-function given by Eq. (151j) leads to a confining theory if 62 > 3/2. 

3 This interval for 62 will be further reduced after imposing the requirements of confinement and the 
infrared singularity being repulsive to physical modes, cf. Appendix B. 



(3(a) = -b 2 a + 



b 2 a + 



--6oK + 
a 



2a 



^,2 



a 



61 a 3 
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To see that we need only to study the IR behavior of the corresponding function X(a) 
cf. Eq. 017]). In this limit 

3 



lim I X(a) H — ) log a = lim ( — — H — ) logn < 

a— s>oo \ 2 I a— >oo \ 3 2 



&9 > - 



2' 



(53) 



which constitutes the general criterion for confinement in its version of the /3-function, 
Ref. [7J. 

The dilaton potential corresponding to the parametrization of /3(a) is 
V(a) = V ( 1 



/3(a) 



3a 



86 2 

a\ — 



a 



((27 - 3)b 2 + 46 a + 2M 2 ) 



■Exp 
■Exp 



•Exp ( — - ■ ExpIntegralE 
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-e"^ I 3b 2 - 46 a - 26 x a 2 + (— - 2b - 2&ia)a 
9 V a 

86 2 









) 


L a. 





(54) 



In this expression, 7 is the Euler's constant, and ExpIntegralE (n, z) = J-j 00 ^-dt is the 
exponential integral function, 60 and 61 are the coefficients appearing in the QCD /3-function 
given by Eq. fl39|) . while Vo, a and 62 are undetermined constants. An interesting exercise 
is to expand Eq. f|54|) in the UV. The result is 



V(a) = ^ |l + |6oa+^-(23^ + 366i) 



a 



2187a 3 
+0(a 4 



[546 2 - 1626 a - 3246!a 2 + 206^a 3 + 1896 &ia 3 ] 



a 



(55) 



The leading orders are determined by the U V parameters bo and b\ , and the unknown con- 
stants 62 and a start contributing at 0(a 3 ). The same feature is shared by the parametriza- 
tion of Ref. [15], although in this reference the order a 3 is replaced by a 8 / 3 . The behavior 
of the dilaton potential for larger values of a can be phenomenologically determined by 
fitting to the heavy quark potential which we will do in the next section. 



5 Fit of Parameters to QQ Potential and Running Coupling 

With the modified dilaton potential of Eq. ( 1541) we can calculate the heavy quark-antiquark 
potential. In Ref. [1] the heavy QQ-potential was in very good agreement with the Cornell 
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5 Fit of Parameters to QQ Potential and Running Coupling 



potential. Now, we recalculate it with the modified dilaton potential V(4>) to focus on two 
particular questions: Firstly, how does the potential from string/gravity theory compare 
with the potential from three loop perturbation theory? This serves as a test whether 
the improvements on the /3-function pay off in the UV behavior of observables. Secondly, 
can the long distance string tension help us fix the remaining parameters? In Ref. [1] the 
calculation was done in the bulk z-coordinate, but in this section it is better to work with 
the variable a = e^. A general derivation for the heavy QQ-potential using a = has 
been given in Ref. [33]. In the following we will add the correct short distance and long 
distance analysis for the first time. We refer to Appendix C for further discussions and 
computation in the small separation limit. 

The first step is to derive the explicit form of the metric which is consistent with the 
/3-function given by Eq. f )5T|) . We work in coordinates dependent on the running coupling 
a as a variable, instead of z or u. From the domain wall coordinates relation e A dz = du, 
one may easily derive 



I - u*^ — - . 

acj) 

Then using Eqs. (jUJ), (jUj) and (gHJ) and solving Eq. we get 




(56) 



da 

dz 




(57) 



with £ given by 



6 



(58) 



while A and e 



are functions of a given by 




(59) 



and 




(60) 



with the fixed constants a* and A* defined at the energy E: 



a, = 0.25241, E = e A *A = 3 GeV. (61) 

Then the new metric in Euclidean space, which includes the new warp factor h(z), follows 
from Eq. <^ using Eqs. (jT5l) and (|57|) . and it reads 
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ds 2 = h(z(a)) 



1 



-(-d1? + d&) + e% -Pe^da 2 



(62) 



(Az(a)) 



with h(z) given by 



h(z) 



(63) 



One advantage of the present computation starting from the /3-function given by Eq. ( ]5Tj) 
is that we need a* as single input value for a, in contrast to the two values we needed 
within the formalism based on the warp factor h(z) of Sees. |2]and[3j cf. Eq. (|29|) . § 

Note that in A(a) the /3-function enters the integral in the denominator, and in e D (a) 
it appears once in the denominator and once in the numerator of the exponential function. 
Taking into account that (3{a) < one finds that both A and e D become large for small a. 
For large a, e D becomes small, while A behaves as — jMoga in this regime. 

The general procedure to compute Vqq within the classical approximation is similar 
to the one used in Refs. [H [13]. The heavy quark potential follows from the Nambu-Goto 
action for a rectangular Wilson loop with a short spatial side and a much longer time side, 
i.e. 



The picture is given by a string stretched between a quark and an antiquark, located at 
x\ = | and X2 = — f respectively, which dips into the bulk of the background AdS^-space. 
The separation p between two quarks as well as the potential energy can be expressed 
as functions of cto, which is the value of a at the mid-point between the quark and the 
antiquark. 

In principle the Nambu-Goto action Sng can now contain a new string length l s com- 
pared with the Nambu-Goto action defined in section [TJ 

4 With the /3-function of the previous model we cannot define e D in the same way because we effectively 
need a cutoff at small a, cf. Eqs. and (|119[) . In this case the computation of da/dz should be done in 
a different way. Using Eqs. (TJU), (TJU) and (j5o) one has 



(W) ~ e 



-T-V 



~ e 



— Sng 



(66) 



da 
dz 



p(a)W(a)e A ^, 



(64) 



where A(a) is given by Eq. (|59|) and 



W{a) = W(a ct )e 




'.a 



(65) 



The cutoff a ct introduces a new integration constant, which is multiplicative and related to a factor e KUV 
in a, cf. Eq. The correct QCD /3-function in the UV, however, renders the integration in Eq. (|65|) 

finite. No multiplicative constant is needed. Only a single input value a* is sufficient. 
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5 Fit of Parameters to QQ Potential and Running Coupling 



where h a b is the new induced worldsheet metric defined by h a b = G^u^r^ir- 

To obtain the heavy QQ-potential, we need to express the separation p between the 
quark and the antiquark, as well as the value of the potential Vqq as functions of ckq, 
which is the value of a at the mid-point between the quark and the antiquark. Similarly, 
A = A(oto) and O = 0(oiq) are the values of A and at the mid-point. The separation 
between the quark and the antiquark is given by [13J 



p{a ) = 2£e~ Ao • / - f g ~ =da , (6£ 

with 



o Vl - a~^e~ 4A 



A = A-Ao, (69) 
= 0-0o, (70) 
a = — . (71) 

The bare potential Vqq calculated from the Nambu-Goto action is divergent, so we have 
to regularize it. The divergence means that the quark- antiquark pair becomes infinitely 
heavy. An obvious way to remove this divergence is to subtract the rest mass of the two 
heavy quarks [T3J [16] . Then the finite part of the potential is 



_ 4 

V QQ {ao) = 



_ 4 



,Y3e D+A (l- Vl- a-le- 4A ) r oc 

da V • L - daal- e D+A 

o vl-a"le- 4 ^ Jao 



(72) 



Combining Eq. ( 168|) with Eq. ( I72|) . we obtain the heavy quark potential as a function of 
the separation between the quark and the antiquark. In the numerical computation of the 
second integral of Eq. (172|) we replace the variable a — > so that the integral transforms 
into an integration between and 1 for the variable a, which is much easier to compute. 
The regularization procedure ensures that the integrals are ultraviolet convergent as can 
be easily proved. 

In order to fix the three parameters Vq, a and b 2 of the dilaton potential and the string 
constant l s , we will study separately the short distance, i.e. the ultraviolet (UV) regime, 
and the large distance, i.e. the infrared (IR) regime. The parameter Vq is relevant in the 
UV, while a and b 2 become important in the IR, and l s naturally scales the potential. 
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Let us first study the infrared properties of integrals appearing in Eqs. (168]) and ( 172]) . 
For this purpose we focus on the large a behavior of the /3-function, where —b 2 a is the 
relevant term, cf. Eq. (15TT) . As we show in Sec. H] a value b 2 > 3/2 ensures that the theory 
is confining. On the other hand, in the limit a — > 00 the argument inside the square root 
in Eqs. f )68|) and (172]) becomes 

l-a - §e -4A ~l-a^, 0<5<1, (73) 

which is negative for b 2 > 3/2. So, a value of b 2 > 3/2 means that «o cannot exceed some 
upper limit «q, and in the limit a — > a$ then p diverges. As an example, a$ = 1.38 for 
b 2 = 1.7, and a^ = 1.08 for b 2 = 2.3. Note that these upper limits are not too high, taking 
into account that the PDG data approximately relate a = 1.45^9 43 to E = 0.6 GeV (see 
Eq. fl52|) and discussion). 

The condition that the second integral in Eq. ( 1721) is finite upon integration to infinity 

f°° 1 -L-i 114 

/ da — a b 2 3 a"^~9 62 < 00, (74) 

Ja b 2 

necessitates | — 5- — |&2 < —1, which is fulfilled only if 62 7^ 3/2 and positive. The 
integration (1741) is convergent, but this convergence is very slow for the values of b 2 we 
consider here. As an example, the convergence for the value b 2 = 2.3 is reached only when 
a ~ 10 20 ao- In practice we handle this problem by analytically computing the integral for 
large values of a, where (3(a) ~ —b 2 a. This approach is excellent for the interval (5a , 00). 
For values close to ao, i.e. in the interval (ao,5«o), we perform a numerical integration. 

The parameters a, b 2 and l s must be chosen to reproduce the physical value of the string 
tension a = (0.425 GeV) 2 . From a numerical computation of the heavy QQ-potential in 
the regime p ~ 5 GeV -1 , we find that these three parameters are constrained according to 
the relation 

— = 3.51GeV-[ s , (75) 
a 

as one sees in Fig. [71 

The parameter l s then follows from Eqs. (1521) and ( 175|) . and we get: 

l s = 1.45 GeV -1 . (76) 

This value is different from the string length Z s = 2.62 GeV -1 used in Section [TJ with the 
guessed metric. Clearly a readjustment of the form of the metric may also lead to a read- 
justment of the string length. At this point we argue that the value of l s is unambiguously 
fixed, and it is not possible to accommodate a value of T s equal to l s within our present 
analysis. 
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5 Fit of Parameters to QQ Potential and Running Coupling 




Fig. 7: Functional dependence of b 2 /a versus l s constrained by the condition that they give a 
physical value for the string tension a = (0.425 GeV) 2 . We consider in this plot the region of the 
parameters set close to b 2 = 2.3. The line corresponds to Eq. (75\ ). 

From our previous analysis and from Appendix B we see that a value 1.5 < 62 < 2.37 
and 0.29 < a < 0.47 satisfying Eq. f )52|) ensures that the theory is confining and the running 
coupling is well reproduced. Even when these intervals are very narrow, one can desire to 
get concrete values for b 2 and a. To this end we study the lowest ++ and 2 ++ glueballs. 
In the presence of both a gravity field and a dilaton field a careful separation of the scalar 
degrees of freedom has to be made [T7J [18]. We take from the second reference [18] the 
corresponding effective Schrodinger potential, which is given by 

V t s ^(z) = (B' t (z)y + B';(z), i = 0,2, (77) 
where the functions B (z) and B 2 (z) differ for the ++ and 2 ++ glueballs 



Bo(z) = \A{z)+ l -\og[X\z)}, (78) 

B 2 (z) = ~A(z). (79) 

A(z) is the Einstein frame scale factor and X[z] = X[a(z)] has been defined in Eq. f )47|) . 
being the dependence a(z) given by Eq. fl57|) . Then we can solve the Schrodinger equation 

^ n {z) = m 2 J, n {z) . (80) 

Best values for the glueballs m ++ = 0.921 GeV, m 2 ++ = 1.462 GeV come out too low for 
62 = 2.3 and a = 0.45, where b 2 and a are constrained themselves according to Eq. ( 152]) . 



o 2 



dz- 



+ V 



Schr. 



» 
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These values for the parameters, which are close to the limit of good infrared singularity, 
cf. Appendix B, give optimum results for the glueball spectrum. The same happens in the 
fit of the running coupling in Sec. HI This justifies the choice 

62 = 2.3, (81) 
a = 0.45. (82) 

It is interesting to note that the glueball spectrum probes larges values of a in the dilaton 
potential than what the heavy QQ potential does. For instance, at the energy of the ground 
state for ++ , the potential is tested up to values of a = 11.1, and for 2 ++ up to a = 16.4. 

The resulting running coupling follows from the /3-function of Eq. (IBTjl by just using 
a single input value (loTT) . See the discussion after Eq. ( 1631) . The behavior of the running 
coupling and its comparison to data from PDG is shown in Fig. [HJ The PDG values for the 
strong coupling are obtained via the PDG web tool [8]. To check the numerical consistency, 
we compare in Tab. [1] the values of the running coupling at several energy scales before 
and after the modification of the dilaton potential, with the corresponding values from 
PDG. The AdS/QCD model can match perturbative QCD-calculations with a very good 
accuracy of about 1%. This is important if one wants to connect a perturbative Monte 
Carlo cascade with non-perturbative QCD physics in parton fragmentation. 



Energy Scale 


a s 


+ 




a 


a 


[GeV] 


(PDG) 






Model 1: h{z) 


Model 2: h(z) 


12 


0.16907 


0.00436 


0.00426 


0.16179 


0.16809 


8.2 


0.18463 


0.00525 


0.00510 


0.18431 


0.18492 


5 


0.20994 


0.00688 


0.00662 


0.21371 


0.21279 


2.6 


0.26708 


0.01178 


0.01106 


0.26809 


0.26640 


2 


0.29942 


0.01522 


0.01408 


0.31023 


0.29668 


1 


0.4996 


0.05799 


0.04598 


0.82393 


0.42618 



Tab. 1 : Values of the running coupling at different energy scales compared with PDG data fE^. 

The values of the running coupling provided by the old warp factor h(z) of Eq. (|2]) 
are doing quite well in comparison with the ones given by PDG within the errors, but the 
values obtained from the new dilaton potential (|54"|) are closer to the experimental data. 

Then only one parameter Vq or I = ^=gpF remains to be fixed. Focusing on the UV 
regime, we can perform an analytical study of the short distance regime by expanding 
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5 Fit of Parameters to QQ Potential and Running Coupling 




Fig. 8'. The running coupling as a function of energy scale. We show as a full (blue) line our 
result stemming from the parametrization of Eq. \51\) , with the values of the parameters given by 
Eqs. \81}) and ItSfy) . The interpolation of the PDG data is displayed as dashed (black) lines, with 
its corresponding error. The point shows the input condition, Eq. h61\) . 



Eqs. (|68|) and ([72]) in powers of oto- The details of the computation are provided in Ap- 
pendix C. The result in NNLO is 



2Pai'\p) 



v qq{p) 



0.359 + 0.5336 a (p) + (1-347^ + 0.692&i)ajj(p) + O(a s ) \ . (83) 



QQKrj nil P 

By reversing Eq. f )68|) . cf. Eq. (I152p in Appendix C, we get the following functional form 
for cto as a function of the separation p 

1 



a (p) 



0.0796^ + 



b log 



1.32£ 



O log- 



+ ^log(6 log M ' 32 ' 



1.32£ 



b log 



1.321 



• £log(Mog' 



P 



(84) 



The heavy quark potential can be directly compared with perturbation theory (PT). 
Vqq is computed in PT as an expansion in powers of the QCD running coupling [19]. It 
has the form 

N 2 c -l a v (p) 



Vpt(p) 



a v = a PT < 1 + (oi + 47T7£& 



2N r 



P 



5) 



where up to the third order 



47T 



+ 



47T7 S (2ai6 + 47r6i) + An' 



7T 



- He )K + «2 



16tt 2 



*6) 
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In this expression b and bi were defined in Eq. (139|) . a PT is the perturbative QCD running 
coupling, and the coefficients a% and a 2 were calculated by Fischler [20] and by Peter [21] 
and Schroder [22], respectively. We use the convention of Ref. [19] and Ref. [22J: 

«i = ^N c - ^N f , (87) 



/4343 , 2 vr 4 22 A Ar2 /5081 . 16 



a 2 = + ^ + — C(3) )Nt - + — C(3) iV c iV> 

2 I 162 4 3 V 1 6 V 324 3 SV ; ' f 



/55 A „ n .\N f 100 Ar2 
+ 4((3) — i H Nf. 



Practically, we use the two- loop running coupling constant, which has the following form 

a PT {p,d) = j—r , (89) 

& o-log(^)+fe-log(21og(^)) 

where d is an undetermined parameter which relates the scale \x of the running coupling 
with the distance p, i.e. p = d/p. 

The expansion ( 183]) is similar to that of PT, Eqs. fl85|) -f l86|) . except that there is an extra 

1/3 

power a at every order. This is a common prediction of all the renormalization group 
revised models constructed by the general procedure of Kiritsis et al., cf. Refs. [HE]. At 
first sight this difference is a matter of concern. In order to fix the unknown parameter Vq 
of the dilaton potential, we proceed in the following way: We find that Eq. (153"]) depends 
on the factor £ 2 /l^ = —12/(Vq1 2 ), and a numerical comparison between the leading orders 
in Eqs. ( l8~3~]) and ( 185]) in the regime 0.06 GeV" 1 < p < 0.20 GeV" 1 gives us the value 
— VqI 2 = 1.31, which is confirmed in a big range of l s . Using the determined value of l s we 
get: 

V = -0.623 GeV 2 . (90) 

A convenient choice for the parameter d follows from a direct comparison between the 
argument inside the logarithms in Eqs. ( I84p and 



d = 1.32M = 1.53. (91) 

We choose the value of A given in Eq. (]3]). Another possibility would be to use a A which 
follows from QCD studies for the running coupling. In either case the value of d is chosen 
in such a way that it compensates a change in A. 

We compare in Fig. |H]the numerical result for the heavy quark-antiquark potential from 
Eq. (172]) . with the perturbative result of QCD given by Eqs. ( 185]) and ( 186]) . in the small 
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distance regime 0.06 GeV -1 < p < 0.20 GeV" 1 . The upper bound in p is motivated by 
the fact that our result for the running coupling fits well the experimental data in the 
regime E > 4 — 5 GeV. We also show in dashed lines the short distance expansion of 
Eq. (1831) up to leading order 0(a^ 3 ), next-to-leading order 0(a^ 3 ) and next-to-next-to- 
leading order 0(al°^ 3 ). We see that the comparison with the numerical computation of 
Vqq and perturbation theory is quite accurate, although it seems to be that the series f )83|) 
is slowly convergent. It is important to note that the full string potential is rather close 
to the perturbative potential in spite of the expansion containing different powers in the 
Yang-Mills coupling. 



\p-V QQ (p)\ 
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Fig. 9: The heavy quark- antiquark potential as a function of the distance p for small p. The result 
stemming from the dilaton potential of Eq. |5^P is shown as a full (blue) line. It follows from a 
numerical computation of Eqs. [68\) and (Tip. The perturbative computation, Eqs. {85\) and [86}) . 
is displayed as a dashed (red) line. The short distance expansion, Eq. (S3) , is displayed up to 
leading order (TO), next-to-LO and next-to-next-to-LO, as dashed (green) lines from bottom to 
top, respectively. We consider in this plot l s = 1.45 GeV^ 1 and Vq = —0.623 GeV 2 . 



A full numerical computation gives the heavy QQ-potential from Eqs. f[68|) and ( 172]) . 
The result is shown in Fig. [TU] . One can try to fit our result with the Cornell form of the 
potential: 

v§r 11 (p) = -- p +°-p+ c > ( 92 ) 

and obtains the values 

a = 0.42, a = (0.415 GeV) 2 , C = -0.14 GeV, 

which are rather close to the accepted values [231 El] , arid also to the values obtained in 
Ref. [Tj. In the numerical computation of Eq. f l72|) we have conveniently normalized the 
result by adding a constant C in order that the QQ-potential vanishes close to p = 2 GeV -1 . 
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Fig. 1 0: The heavy quark- antiquark potential as a function of p. Our result stemming from the 
dilaton potential of Eq. |5^| j is shown as a full (blue) line. It follows from a numerical computation 
of Eqs. (ESP and (7^ ). The dashed (red) line corresponds to the Cornell potential, Eq. (fM). 

6 Dilaton Potential and New Warp Factor 

With the three parameters Vo, a and 62 determined in the previous section we can obtain 
the dilaton potential Eq. ( 154"|) . the running coupling and the modified warp factor h(z) of 
Eq. (153]). 

The parameters governing the UV-asymptotic behavior of the dilaton potential are 
bg = |^ and &i = -^2, the known coefficients of the perturbative /3-function, in addition 
to Vo, as can be seen in Eq. ( |55l) . So, the potential is consistent with the QCD /3-function 
at two-loop level. 

The dilaton potential is plotted in Fig. CD] and compared with the one obtained from 
the warp factor h(z) of Eq. (T5]). The comparison shows that the modified dilaton potential 
becomes slightly flatter. This flattening is in line with the result given by Ref. |15j . although 
the parameters provided in this reference together with our value of l s produce a potential 
of the order of 10 4 GeV 2 in the regime of physical interest, a ~ 0.3, in contrast to the value 
~ 1 GeV 2 given by our potential. 

It is interesting to compute the warp factor h(z) that follows from the new dilaton 
potential of Eq. (|54"|) . This way we can close the circle of investigation in the present 
paper. The mathematical procedure is as follows: From Eqs. floTl) we get a first order 
differential equation 

da = V-3V cA(a) _ D(a) ^ ^ 
dz 6 

where A (a) and D(a) are given by Eqs. ( 159]) and ( 160|) respectively. This equation can be 
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Fig. 1 1 : The dilaton potential as a function of the running coupling a. The result corresponding to 
the parametrization of Eq. 155\) is displayed as a continuous (blue) line. We show for comparison 
as a dashed (red) line the potential stemming from the warp factor h(z) of Eq. computed in 
Sec.\M 

solved by imposing the boundary condition a(z — > 0) — » 0, due to the asymptotic freedom. 
The solution is delicate, because there are problems to find a numerical solution of Eq. ( 193]) 
near the boundary z — 0. At this point the l.h.s. of Eq. ( 193]) is divergent. To overcome 
this difficulty, we proceed in three steps. First, we consider the lowest order perturbative 
expansion of the /3-function and rewrite Eq. fl93|) in the deep UV, 



dz ~ u 6 ^ x 
The solution of this equation is: 



o a^e'D^v . (94) 



a vv {z) = - 1 , (95) 
b log (Az) 



with 



A = = 237 MeV . (96) 

6 

Note that the value of A is rather close to the one determined using h, cf. Eq. ()3]). This 
small discrepancy could be improved when higher orders in the perturbative expansion are 
considered, as it has been observed in the analytical computation of the QQ-potential, see 
Appendix B, cf. Eqs. (fl63l) and ( 1T64j) . 

Choosing z = 1.2 • 10~ 4 GeV -1 we obtain a^o) = 0.0718 as initial value in the deep 
UV to find the numerical solution of Eq. ( 193]) . In the second step we consider three 
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orders in the UV expansion of the ^-function, Eq. ( )5T|) . and A(a), and solve Eq. (1931) from 
z = 1.2 • lO^GeV" 1 to Z! = 0.125 GeV" 1 with 



A(a) = C A + ^~ + % log a + — (&0&2 - 3b 2 a - Qb ha 2 - Qbja 3 ) a + 0(a 2 ) , (97) 
b a Oq oOqCT v 

where Ca = —0.1057 is a constant to make A(a) consistent with the input condition (|6"T1) . 
Finally from Z\ to infinity we use the input functions to solve the equation fully numerically. 
We have checked the stability of the solution changing z and Z\, and considering higher 
orders in the expansion, Eq. (J9~7|) . Once we know a(z), the corresponding warp factor is 
easily computed, and it reads 

h( z ) =e 2A{a{z) \a{z)fz{kz) 2 . (98) 

The strength in the Nambu-Goto action, Eq. (jSJ), is determined by the factor h(z)/l 2 , 
and so it is more relevant to consider this quantity when comparing different models. 
The result h(z)/l 2 is shown in Fig. [12] and compared to the warp factor of Eq. (j2J) that 
was first proposed in Ref. pQ. I 2 is approximately a factor 3 larger than I 2 , and this is 
reflected also in the values of h(z) compared to h(z). The numerical agreement is rather 
good up to z ~ 2GeV _1 , in spite of the fact that h(z) vanishes in the UV. h(z) has a 
singularity at z m = 3.65 GeV -1 , which occurs at slightly larger values than in h(z), for 
which z m = VT^~e/X = 2.73 GeV" 1 . 

In the deep UV, i.e. at leading order in the perturbative expansion of the /3-function, 
the analytical solution of Eq. (193]) has a simple form: 



A\ 2 i 



M°0 = \j) (99) 



A x 2 
A 



(-60 log (Az)) 



(100) 



This functional form of the calculated warp factor is very similar to the guessed warped 
factor besides the power 4/3. 



7 Discussion and final remarks 

The analogy of the bulk coordinate z with the inverse energy resolution has triggered the 
guessed warp factor in Refs. [H [2] , which was based on a naive equivalence with the running 
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7 Discussion and final remarks 




Fig. 12: The warp factor divided by 1% as a function of z. The factor of Eq. p|) proposed in 
Ref- W is shown as a dashed (red) line. The full line corresponds to the warp factor that follows 
from the dilaton potential of Eq. |5^| j using the parameters determined in Sec. 0. 

coupling of QCD. A careful analysis of the resulting dilaton potential gives the evolution of 
the dilaton field in the bulk and consequently the running of the QCD gauge coupling. The 
infrared physics of the model of Refs. [U [2] was satisfactory to fit the string tension, but it 
fails to give a good UV-behavior for the /3-function. Because of the second order Einstein 
equations and the correlated behavior of the dilaton in the infrared and ultraviolet, which 
is not constrained by the correct QCD /3-function, one needs two boundary conditions 
to interpolate the gauge coupling between the charmonium and bottonium masses. This 
feature weakens the idea of holography for the old ansatz, which determines the field 
theoretic behavior of our 4-dimensional world from the physics of gravity in 5-dimensional 
anti-de Sitter space. 

In order to have the Yang-Mills theory as a 'hologram' of the physics happening in five 
dimensions we assumed a new ansatz that improves the UV-behavior by using the QCD 
/3-function as a constraint. Thereby we found a dilaton potential which is consistent with 
QCD in the UV region. The resulting short distance heavy quark potential r ■ Vnn{r) has 
a similar shape as the 3- loop expression derived by Brambilla et al. [19J. The numerical 
comparison of the QCD and string heavy QQ-potential is rather good in spite of the fact 
that the leading term in string theory proportional to a 4//3 deviates from the QCD-potential 
proportional to a. By calculating the NNLO-expansion, we show that the expansion in 
a is slowly converging. With a new string length l s in the Nambu-Goto action we can 
also match the long range QQ-interaction. We have closed the circle of investigation by 
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computing the warp factor corresponding to the new ansatz, and found that the scaled warp 
factor is similar to that of Refs. [TJ |2] in the region of interest 0.5 GeV^ 1 < z < 2 GeV -1 . 

In the procedure we have proposed, we fix three parameters (62, ol, Vq) corresponding 
to the dilaton potential and one parameter l s equal to the string length by using three 
constraints, namely a good behavior of the QQ"P°t;ential in the IR and in the UV, and 
a good behavior of the running coupling in the regime 0.6 GeV < E < 15 GeV. New 
parameters may have to be included in order to describe the glueball spectrum or further 
observables with precision. One important point of this analysis is the general criterion for 
confinement of Ref. [7] and the requirement that the infrared singularity being repulsive to 
physical modes, which help us to set a narrow window in our parameter set, in particular 
for the parameters controlling the infrared behavior of the theory (6 2 , «)• 

Recently the question "How Well Does AdS/QCD Describe QCD?" [25] has been asked, 
and depending on the feature the answer varied - the accuracy was estimated between 10% 
and 25%. In our case of pure gluon dynamics, we have shown that the accuracy is much 
better, and therefore we look optimistically towards further tests of the action in finite 
temperature calculations [26] . 
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Appendix A: Infrared and Ultraviolet Properties of the Gravity Dual 
Theory 

In this appendix we study the infrared (IR) and the ultraviolet (UV) properties of the 
5-dim Nambu-Goto theory with the "guessed" metric h(z)/(Az) 2 , cf. Eq. (J2]), which we 
develop in Sees. |2] and El We apply technical details of Refs. |U [7], which are shortly 
reviewed in Sec. HJ The warp factor h(z) has a singularity at z m = a/1 — e/A, so the bulk 
coordinate z is restricted to z < z m . The IR expansion of h(z) is 

log(i) 1 (l-2e)log(i) 

w = *m<* -*T4 +ciQt {_M) ' (101) 

with 

i = z m -z. (102) 
Using Eqs. ( JT5]) and f[2"Tj) we get a second order differential equation for </>(£), which writes 



(103) 



This equation can be solved for several orders in the IR expansion of h(£). The theory 
of differential equations gives the general solution by adding to the special solution of the 
inhomogeneous equation the full set of homogeneous solutions. The result is 

<KO = ^(logO 2 + ci log£ + c 2 + 0(£log£) . (104) 
lb 

The parameters Ci and c 2 are two unknown constants corresponding to the homogeneous 
solutions which have to be fixed by two conditions. Obviously the three terms that we 
show explicitly in Eq. ( 1104jl correspond to the lowest orders of the homogeneous and 
inhomogeneous solutions. At this point it is preferable to write the solution in this way 

^ = Ta ( lo S —) + + lo § > ( 105 ) 
16 V w m J 

where the constants u m and related to c\ and c 2 as 

8 4 
w, R = e"3 Cl , «m = --cf + c 2 . (106) 

Setting w IR corresponds to setting the scale. We can solve Eq. f j2Tl) numerically for the 
full range of z < z m using the IR behavior of Eq. fl 105ft as a boundary condition. In the 
concrete calculation we choose 0(£i) and 0'(£i) whh very small. Doing that, we have 
checked that the result of Fig. [Tj is exactly reproduced for 

co m = 4.55 GeV -1 , n m = -0.758 . (107) 
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These numbers are stable in the deep infrared near the singularity, £1 ~ 10 7 — 10 4 GeV . 
So, this choice of the constants u m and K m is equivalent to the boundary conditions 
of Eq. (EE]). 

The IR behavior of A(^) can be obtained from Eqs. f|T5|) . (11 01 1) and (11051) . and it reads 

4« - -s (* 2 -\ £) (io^fcx) + °« logf » • ^ • 

(108) 

Once that we know A(£), the IR asymptotics of the superpotential W can be computed 
using Eq. (jUJ) and taking into account the domain wall coordinates relation e A dz = du. 
It reads, 

W(a) = W OQ al (loga^e^^ ( 1 - 3 + ^ l — + O ((loga)- 1 ) | , cw oo , 

\ 2v/3 (log a) 3 U 7 

(109) 

where the constant is 



3 6(1- e)2 A , 

4y tog (i) u m 

The dependence in k ir appears at 0((loga) - ^) in the bracket of Eq. (I109p . The fact that 
the superpotential W(a) grows faster than «3 in the infrared ensures that the theory is 
confining and also that there is a mass gap in the spectrum [TJ [9] . Note that this is true 
independently of the values of u m and K m . 

The dilaton potential can be obtained from W using the third relation of Eq. (146p . or 
equivalently from A(z) and Eq. fl20|) . In this regime it behaves as 



V(a) = VooctQoga) ( 1 - 2(2 + Km) - 1 + O ((logo:)- 1 ) 

V Vo (log a) 2 



a — > oo , 
(111) 



where the constant Voo is 

16 r, . 

Ko = -y^. (112) 
The general form of the potential in the infrared that has been studied in Refs. [U [7] is 

V(ce) ~ a 2Q (loga) p , a-^oo. (113) 

The solution of Eq. f illip corresponds to Q = 2/3 and P = 1 in this notation, in addition 

4 /j 

to an extra factor oga which is dominant with respect to log a, but subdominant with 
respect to ai This subdominance, in addition to the fact that Q < 2y / 2/3 ensures that 
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the IR singularity is of the good kind according to the criterion of Giirsoy et al. |7J , which 
means that the singularity should be repulsive to physical fluctuations. To avoid any doubt 
about the extra factor, it is possible to prove that it can be removed in Eqs. (11091) and 
( lllip by shifting the dilaton field 

^=0U^, (114) 

which doesn't affect the IR assymptotics. The resulting expression for the assymptotics of 
the dilaton potential is: 

4 

V(a) ~ «s log a , a — > oo . (115) 

Note that the asymptotics of Eq. (11081) has not been studied by the authors of Ref. [UH], 
and so this case is not listed in Table 1 of Ref. [7] . One can add our result to this table. 
The IR behavior of X(a) can be computed from Eqs. (1451) and (I109|) : 

Yf s 1 V3 1 3 1 y/3(3 + 2k ir ) 1 

X(a) = ; : =- + ..., a— too. 116 

2 4 (loga)5 8 log a 16 (loga)l 

The term oc l/(loga)^ comes from the extra factor e^^" in W(a). In our case, the 
limit 

lim (X(a) + - J log a = K, (117) 



a— >oo V 2 

leads to K = — oo. 

One important aspect of this analysis is that one has to fix two integration constants, 
u m and K m , using initial conditions. This contrasts with the analysis of Refs. [3 |9], where 
they show that just one initial condition is enough. One possibility studied in Refs [TJ [9] 
is fixing one of the parameters by requiring that the bulk singularity is not of the "bad 
kind", which means that the singularity should be repulsive to physical fluctuations, cf. 
Eq. (E.28) of Ref. [9]. In our case the singularity of our solution at z — zir is of the good 
kind, independently of the values of u m and k ir , cf. Eqs. (I109p . (lllip . (I112p and compare 
with Eqs. (E.27) and (E.29) of Ref. [9|. To give a new perspective to this issue, we can 
study the UV behavior of the dilaton using h(z) of Eq. ([2]). Following the same procedure 
that we explained for the IR but considering an expansion at small z, one gets 

0(*) = -^ + k uv + %4tT— + ^ 2 )' *">°> ( 118 ) 
z elog(i)w uv 

where u vv is a parameter setting the scale and k vv is another parameter. They play the 
role of u m and k ir respectively in the UV. In fact these parameters are related, in the sense 
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that setting k uv in the UV then sets n m in the IR (the same for co m and u vv ), cf. Sec. [3j 
From Eq. (11181) and using the same procedure that we explain above, one gets 

at \ 3 9 a 9(3 + 2k uv ) a ( a \ 

3(a) = --a - — + O — , (119) 

v ' 2 41oga 8 (loga) 2 \(\oga) 3 J 

which has the drawback not to be consistent with asymptotic freedom as found in QCD 
and has motivated our improvement in Sec. HI Note that the parameter k vv enters in 
the expansion of (3(a), Eq. f l 1 1 9 j) . A well defined function (3(a) from perturbation theory 
doesn't have this parameter, and this is the starting point of the program followed in 
Refs. [U [7J [9]. As we explain in Sec. [3] we fix /% v and the scale parameter w uv (equivalently 
k 1k and cu IR ) by using two input values for a, cf. Eq. 



Appendix B: Infrared Properties of the Improved Gravity Dual Theory 

In this appendix we study the infrared properties of the 5-dim Nambu-Goto theory 
with the improved metric h(z) / (Az) 2 , cf. Eq. ( 15^]) . which is proposed in Sec. [Hand further 
developed in Sees. Eland El The ultraviolet properties are fully dictated by the UV behavior 
of the /3-function, so we will just focus on the IR asymptotics. As we will see later, the warp 
factor h(z) has a singularity at some finite value z lR , so that the coordinate z is restricted 
to z < z m , like in the guessed metric h(z)/(Az) 2 studied in Appendix A. 

The IR expansion of A(z) can be computed from the function A(a) and the functional 
dependence a(z). In the IR, i.e. a — > oo, Eq. (T59]) can be integrated out using the 
expression of the (3- function given in Eq. (J5TJ) . One gets 

A(a) = A(a*)-^-\og(—)+... a oo . (120) 

The computation of e D involves an integration from to a, and so it must be performed 
more carefully in order to retain the UV convergence, cf. Eq. (1601) . Inserting the full 
expression of (3(a) into the integrand of Eq. ( 1601) . and considering an expansion at large a, 
one gets 



2 a. 

— e « 
9 



(6 2 - 2a(6 + ha)) - + O(a ) 
a 



a ^ oo , (121) 



where for simplicity we have defined the constant 



Q) = 7: [62 (3 - 2 lE ) - 2a(2b + b x a)] 



122) 
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The integration is convergent in the UV, and the dominant contribution in the IR comes 
from the logarithmic term in the r.h.s. of Eq. (I12ip . From Eq. (!60|) . and using Eqs. ( 15T|) 
and (11211) one can easily compute the IR asymptotics of e D . It reads 



Co /a ,-lb 2 



[l + £>(e~§ 



a — > oo . 



(123) 



The functional dependence a (z) can be computed from Eq. (!57|) . which one can write 
in the following way 



dz = £e D - A da. 

Then the function z(a) follows from an integration of Eq. (11241) . 



dz 



a(z) 



D D(a)-A(a) da ^ 



«i = ac(Zi) . 



(124) 



(125) 



Inserting Eqs. (11201) and (11231) into the r.h.s. of Eq. (11251) . we can perform the integration 
analytically. After inversing the solution, one gets 



a{z) 



C*(zi — z) + a x 



b 2 9 



b-2 



b2 



1 



[C t ■ {z m - z)} 



S-b 2 



where we have defined 



b 2 



ih 2 - 1 ' 



and the constant C* is 



z m — z l + 



1 c A(a,)-C a 



Q i 



c* 



z z ra , (126) 



(127) 



'128) 



Inserting Eq. ( 1126}) into Eq. ( 11201) then one finally gets the IR asymptotics of A(z), 
which reads 

A(z) = 5 ■ \og{z m - z) + . . . , z^z m . (129) 
The IR behavior of <f>(z) can be obtained from Eqs. ( l2Ti) and (11291) . and it reads 



--^5(1 + 5) \og(z m - z) + . . . 



z -> z m , (130) 

which corresponds to the asymptotics of log a(z), cf. Eq. (I126p . after taking into account 
Eq. (H2U). 



35 



The procedure to compute the IR asymptotics of the superpotential W and the dilaton 
potential V is explained in Appendix A. The superpotential reads 

W(a) = ^5 -Cl +S -al b2 + ... , a -> oo , (131) 
and the dilaton potential 

V(a) = \cl +2S ■ 5 2 ■ (h\ - 9)a> + . . . , a oo . (132) 

Within the notation of Refs. (HE], cf. Eq. ( 11131) . the solution of Eq. (11321) corresponds 
to Q = |6 2 - As it has been discussed in Ref. [7J, the asymptotics of Eq. (I120p leads to a 
confining theory whenever Q > 2/3, and this is fulfilled in our case for b 2 > 3/2. The IR 
behavior of X(a) reads 

X(a) = —b 2 + ... , a ^ oo (133) 
3 

so that the limit 

lim \X(a) + - J loga = K, (134) 

a— s>oo V 2 / 

leads to if = — oo in the confining case, i.e. b 2 > 3/2, as it was explained in Sec. HJ cf. 
Eq. (J53]). 

Giving the IR asymptotics of Eq. (11131) . a good singularity according to the criterion 
of Ref. [7J is obtained when Q < 2\/2/3. This means that our theory is confining and it 
presents a good singularity for values of b 2 in the range 

- 2 <*2<\. (135) 

Nevertheless, as we will see below the upper bound is too conservative, and a slightly 
larger value is obtained when computing the IR asymptotics of the effective Schrodinger 
potential for the glueball spectrum. To this end, first we will analyze the 2 ++ sector, whose 
potential is 

V 2 s ^iz) = {B' 2 {z)f + B>'{z) , (136) 

where 

B 2 (z) = ^A(z) . (137) 

After inserting the IR asymptotics of A(z), Eq. (11291) . into Eq. (11371) and computing the 
derivatives, one gets the following asymptotics for the effective Schrodinger potential 
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3 „ /3 



^ = 2 , l2 , - 1 j7*r^ + - > 2 ^'" (138) 



The potential diverges to +oo whenever 5 > 2/3, or equivalently 

6 2 <^y|^ 2.37. (139) 

Otherwise the Schrodinger potential would diverge to — oo at z — > z m , so that the 
singularity would be attractive to physical fluctuations. So, the condition Eq. (I139p ensures 
that the singularity is of the good king according to the criterion of Gursoy et al. [TJ. 

The computation in the ++ sector is similar. In this case the Schrodinger potential 
reads 

V (z) = (B^z))* + B£(z) , (140) 

where 



B (z) = ~A(z) + ±log[X*(z)], (141) 



and 



X W = M. (142) 
6a{z) 

Using the asymptotics of A(z), Eq. (I129p . in Eq. (I14ip and the asymptotics of a(z), 
Eq. (I126p . in Eq. (I142p . one gets for V (z) the same IR asymptotics as for V 2 (z), cf. 
Eq. PH]). 

We plot in Fig. [13] the effective Schrodinger potentials for the computation of the ++ 
and 2 ++ glueball spectrum using the guessed value 6 2 = 2.3. 



Appendix C: Heavy QQ-Potential at Short Distance 

In this appendix we give technical details for the analytical computation of the heavy 
quark- antiquark potential in the small separation limit as an expansion in powers of a^. 
The relevant formulas are Eqs. (1681) and (1721) . The result up to NLO has been derived in 
Ref. (EJ. 

In order to compute the series up to NNLO, we need to consider the QCD /3-function 
up to order qjq, so we have to take into account the parameters b and b± (see Eq. fl39|) ). 
We have explicitly checked that the order «q in the /3-function contributes at higher orders 
in the potential. Here and in the following we use the notation a = a ■ a^. The variable a 
is in the interval < a < 1 in the integral of Eq. (|68|) and the first integral of Eq. ( I72|) . 
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Fig. 13: The effective Schrodinger potentials for the ++ (full red line) and 2 ++ (dashed blue 
line) glueball spectrum using hi = 2.3. We also show in full green line the running coupling. 



First we will consider the computation of p{olq) which is given by Eq. ( 1681) . The function 
A(pt) can be computed using Eq. ( 1591) . and the result is 



a 



1 b 

-4(a) = 7— + 4 !og , , , 

b a bg \b + bid 



+ 0(a) 



(143) 



The function e D follows from Eq. ( 16 Up . and after an expansion in powers of «o it reads 

1 



boa' 



(l + diao + ol 2 al + 0{otD) , 



where 



di 

d-2 



4 h' 

J r 77 ) ''ii" • 



9 ' b? 



816g 



(8&£ + 186^6i + 816?) a 2 



(144) 

(145) 
(146) 



The integration in a cannot be computed analytically even after this expansion. In order 
to avoid this problem, we consider the following series 



-3A 



V7 



-4 A 



E 



(2n)(2w- jj|L §(l-nUl-4n)i 



^ (2n- l)(2n)!! 



(147) 



Note that A = A(a) — A(o;o) is also a function of the coupling constant cuo, so it makes 
sense to consider an expansion in ao for the above expression. The result is 



e d-4n M = e (l-4n)^ 5 (l ^n)| ^ + ^ + ^ + Q ^ 



(148) 
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where 

<! = (l-4n)(l-fi)§, (149) 
h = _(i_4 n )(i_ 5 )^[^(i + fi )_ 6l (i_ 4n )(i_ & )]. (150) 

In Eq. (11481) there is a dependence in ao that cannot be expanded due to an essential 
singularity. Combining Eq. (I144p and Eq. (I148p . we get 

e D e d-4n)A = JL e (l-4n)^ ~d-4n)|-2 ^ + ^ + ^ + ( ^ + ^ + ^ + 

o a 

(151) 

Using this expression and taking into account Eq. (1681) and Eq. (I147p . we can compute the 
integration in a analytically for every order in ao and n, and the result is a summation of 
terms involving the incomplete Gamma functions. Making a further expansion in ao one 
gets 

p(a ) = 2le~ Ao [po + Pl a + p 2 a 2 + O{a 3 )] , (152) 



where 



Po 



Pi 



P-2 



tl (2 'f""^ 0-596. 
^ (2n- l)(4n- l)(2n)!! 

16, ^ n(n - l)(2n - „ , , . 

0.0860^ + O.I8O&1 . 



We truncate the summation of the infinite series to ^^oooo ... ; by which the residue con- 
tribution can be neglected to 10 -3 relative precision. 

Next, we will focus on the computation of Vqq, which is given by Eq. (I72p . There are 
two contributions which we call V and V s , corresponding to the first and second integrals 
respectively, so we write Vqq = V — V s . Note that for the first one, the expression is similar 
to that of p(a ), and so we can apply the same procedure that we explained above. In 
order to make the integration analytical we consider the series 

e A (l-Vl- a-ie- 4 ^) °° (2r) _ lV , 
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Using Eqs. (172]) . (1144[) . (I148j) and (I154p we can compute analytically the integration in a 
in the same way as we did for p(«o)- Finally we arrive at this result 

_ 4 

la 3 e A ° 

V ( a o) = Ij 2 K + viao + v 2 a 2 + 0{al)] , (155) 



where 



^ 2n-l 

> 77^ TT^-^TT « 0.398 , 

^ {An- l)(2n)!! 

8, ^ (n- l)(2n- 1)!! , , , 

-&o> V: ^777 — rfr « 0.04236 , 156 

Q ^ {An- l) 2 (2n)!! v ; 



n=l 



«2 



2 f2 1 1 ! ! 

81 E ( 4w -n3 ( 2n)!! + n(4 ° n " 53))6 ° + 9{An ~ 1){8n ~ ^ 

n=l 

« 0.06406^ + 0.131&1 . 

The computation of V s (ao) is a little bit tricky because it does not allow similar expansion 
techniques. First of all it is convenient to change the integration limits from (ao, oo) to 
(0, 1), and so we consider the variable replacement a — > Then we get 
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K = ° / «3 . e D+A rfa = / d t ■ e D+A cfo =: — / a ^f{a, a )da , 

(157) 

where we have defined the function 

f(a,a ) = 4e D+A , (158) 

which involves all the interesting dependence in ao- We can expand this function in powers 
of ao by considering its derivatives 

f{a, a ) = /(a, 0) + /'(a, 0)a + i/"(a, 0)^ + O(ag) , (159) 

where /' stands for derivative with respect to the second argument. This expansion can 
be easily done analytically, but the integration in a must be computed numerically. By 
inserting Eq. (j!59p into Eq. (jl57p . and after making the integration, finally we get 

_ 4 

£a s e A ° r i 
V s {a )= ° l + O.8896oa o + (2.176g + l.ll&i)ag + 0(ag) . (160) 

7T/2 L 'J 

Note that the contribution of the substraction term is larger than Eq. f)155p by approxi- 
mately a factor 3 in the regime which we are interested in. Eq. (I160p reproduces well the 



40 



References 



full numerical computation of Eq. (I157P up to «o — 0.19, with an error of 5%. This «o 
corresponds to distance p ~ 0.20 GeV -1 . 

Combining Eqs. (I155p . (II 56ft and (I160p . we obtain for the QQ-P°tential 

_ 4 

£<y 3 e A ° r i 

Vqo(p) = V - V s = ^— 0.602 + 0.8476 a + (2.1066g + 0.9796i)a^ + C»(a|j) . (161) 

vr/2 L J 

From Eq. H152[) and Eq. f II 6 1 j) we can compute the dimensionless quantity pVqq(p) up to 
order Qq, and we get the final result 

Vqq(p) = -H! " ^ [o.359 + 0.5336 «o(p) + (1.3476g + 0.692^)^) + O(a 3 )] . (162) 
ixlj p l J 

The function ato(p) can be obtained by reversing Eq. (I152p . To do that we make first an 
expansion of Aq in ao- The resulting expression is 



1 i h C + §) «o 

«o(p) = T - \ ~ 7 \ 07 (163) 

6 ° lo § (f J - bo lo S (&o log (f) - I log a (p)^ 




where 



c = 2p b b l/b h^ 1.32£, (164) 



with the parameters po an d Pi defined in Eq. (I153p . Then we can use an iteration method to 
obtain the solution order by order. For example, the first order is given by ^60 log (fjj 
If one substitutes this expression in the r.h.s. of Eq. (11631) . then we get the NLO approxi- 
mation. 

Combining Eq. (I162p and Eq. (I163p we have the <5<Q"P°tential as a function of the 
separation p between the quark and the antiquark in the limit of small separations, i.e. 
p — > 0. This result agrees with Ref. [T3] at leading order. We corrected a mistake in this 
reference at next to leading order for V s . Fully NNLO expressions are provided here. 
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